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Abstract 
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by-light contribution to muon g — 2. 
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1 Introduction 



Here, I discuss the contribution to the muon g( — 2 of a hadronic bubble connected 
to the external static magnetic source through one photon leg and to the muon 
line with another three photon legs. This corresponds to the so-called hadronic 
light-by-light contribution to the muon anomaly = {g — 2)/2. Recent reviews 
are in [1,2]. One of the six possible photon momenta configurations is shown in 
Fig. 1 and its contribution to the vertex — |e| u{p') T^{p — p') u{p) Ap is 
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where p3 — > is the momentum of the photon that couples to the external 
magnetic source, q = Pi + P2 + Pd, and m is the muon mass. The dominant 
contribution to the hadronic four-point function 

n'''^"^(pi,P2,P3) = 
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comes from the three light quark {q = u, d, s) components in the electromagnetic 
current V^{x) = ^Q'j'^q (x) with Q the quark electrical charge matrix. 
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Figure 1: One momenta configuration of the hadronic light-by-light contribution 
to muon g — 2. 

Using gauge-invariance, one can write 

W "'^ Pl,P2,P3) = -P3X ^ ^ (3) 

and therefore one just needs derivatives of the four-point function at pa = 0. The 
contribution to a^^ is 

C=4^W(:^ + rn)[7^7l(^ + m)^^}. (4) 



The four-point function U.f"^°'^{pi,p2,p3) is an extremely difficult object involving 
many scales and no full first principle calculation of it has been reported yet. 
Notice that we need momenta pi and p2 varying from to oo. Unfortunately, 
there is neither a direct connection of aj^^ to a measurable quantity. Two lattice 
groups have started exploratory calculations [3, 4] but the final uncertainty that 
these calculations can reach is not clear yet. 

Attending to a combined large number of colors A^^ of QCD and chiral per- 
turbation theory (CHPT) counting one can distinguish four types of contribu- 
tions [5] . Notice that the CHPT counting is only for organization of the contri- 
butions and refers to the lowest order term contributing in each case. The four 
different types of contributions are: 

• Goldstone boson exchanges contribution are 0{Nc) and start at 0{p^) in 
CHPT. 

• One-meson irreducible vertex contribution and non-Goldstone boson ex- 
changes contribute also at 0{Nc) but start contributing at 0{p^) in CHPT. 

• One-loop of Goldstone bosons contribution are 0{1/Nc) and start at 0{p'^) 
in CHPT. 

• One-loop of non-Goldstone boson contributions which are 0{1/Nc) but start 
contributing at 0{p^) in CHPT. 

Based on the counting above there are two full calculations [6,7] and [8,9]. There 
is also a detailed study of the 7r° exchange contribution [10] putting emphasis in 
obtaining analytical expressions for this part. 

Using operator product expansion (OPE) in QCD, the authors of [11] pointed 
out a new short-distance constraint of the reduced full four-point Green function 

{0\T[V^{p^)V^{p2)VP{-{p,+P2+Ps))]\^{ps)) (5) 

when P3 — > and in the special momenta configuration —p^ ~ — p| >> — (pi -|- 
P2)^ Euclidean and large. See also [12]. This short distance constraint was not 
explicitly imposed in previous calculations. 

2 Leading in the l/Nc Expansion Contribution 

Using effective field theory techniques, the authors of [13] shown that leading 
contribution to aj^' contains a term enhanced by a \og^{n/m) factor where n 
is an ultraviolet scale and the muon mass m provides the infrared scale. This 

leading logarithm is generated by the Goldstone boson exchange contributions 
and is fixed by the Wess-Zumino-Witten (WZW) vertex 71*^77. In the chiral limit 
where quark masses are neglected and at large Nc, the coefficient of this double 



2 



Table 1: Results for the 7r°, rj and rj' exchange contributions. 



Reference 


10^° X 
7r° only vr'', t] and 7]' 


[6,7] 


5.7 


8.3 ± 0.6 


[8,9] 


5.6 


8.5 ± 1.3 


[10] with /i2 = 


5.8 


8.3 ± 1.2 


[10] with /i2 = -10 GeV2 


6.3 




[16] 


6.3 ~ 6.7 




[11] 


7.65 


11.4±1.0 



logarithm is model independent and has been calculated and shown to be positive 
in [13]. All the calculations we discuss here agree with these leading behaviour 
and its coefficient including the sign. A global sign mistake in the 7r° exchange 
in [6,8] was found by [10,13] and confirmed by [7,9] and by others [14,15]. The 
subleading /i-dependcnt terms [13], namely, log(/i/m) and a non-logarithmic term 
are model dependent and calculations of them are implicit in the results 
presented in [6-9, 11]. In particular, k{ij,) contains the large Nc contributions 
from the one-meson irreducible vertex and non-Goldstone boson exchanges. In 
the next section wc review the recent model calculations of the full leading in the 
1/A^c expansion contributions. 

2.1 Model Calculations 

The 71^ exchange contribution was calculated in [6-10, 16] by constructing the 
relevant four-point function in terms of the off-shell T^^'y* {pi)'y* {P2) form factor 
J-'{pI,pI) and the off-shell 7r°7*(g)7(p3 = 0) form factor J-'{q^, 0) modulating each 
a WZW 7r°77 vertex. In all cases several short-distance QCD constraints were 
imposed on these form-factors. In particular, they all have the correct QCD 
short-distance behaviour 

HQ',Q')^^2 ^(g^o)^^ (6) 

when is Euclidean and are in agreement with 7r'^7*7 data. They differ slightly 
in shape due to the different model assumptions (VMD, ENJL, Large NxQM) 
but they produce small numerical differences always compatible within quoted 
uncertainty ~ 1 x 10"^° -see Tabic 1. 

Within the models used in [6-10, 16], to get the full contribution at leading 
in 1/Nc one needs to add the one-meson irreducible vertex contribution and the 
non-Goldstone boson exchanges. In particular, in [8, 9] the one-meson irreducible 
vertex contribution below some scale A was identified with the ENJL quark loop 
contribution while a loop of a heavy quark with mass A was used to mimic the 
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Table 2: Sum of the short- and long-distance quark loop contributions [8] as a 
function of the matching scale A. 



A [GeV] 


0.7 


1.0 


2.0 


4.0 


10^° X 


2.2 


2.0 


1.9 


2.0 



Table 3: Results for the axial- vector exchange contributions from [6,7] and [8,9]. 



Roforoiic-os 


10^-" X 


[6,7] 


0.17 ± 0.10 


[8,9] 


0.25 ± 0.10 



contribution massless QCD quark loop above A. The results are in Table 2 where 
one can see a very nice stability region when A is in the interval [0.7, 4.0] GeV. 
Within the ENJL model, the ENJL quark loop is related trough Ward identities 
to the scalar exchange which we discuss below and both have to be included [8,9]. 
Similar results for the quark loop below A were obtained in [6, 7] though these 
authors didn't discuss the short-distance long-distance matching. 

The exchange of axial- vectors and scalars in nonet symmetry -this symmetry 
is exact in the large limit, was also included in [8,9] while only the axial- vector 
exchange was included in [6, 7] . The result of the scalar exchange obtained in [8] 
is 

a^(Scalar) = -(0.7 ± 0.2) x 10"^° . (7) 

The result of the axial- vector exchanges in [6, 7] and [8, 9] can be found in Table 
3. 

Melnikov and Vainshtein used a model that saturates the hadronic four-point 
function in (2) at leading order (LO) in the 1/Nc expansion with and axial- 
vector exchanges. In that model, the new OPE constraint of the reduced four- 
point function found in [11] forces the vr°7*(g)7(j93 = 0) vertex to be point-like 
rather than including a 0) form factor. There are also OPE constraints 

for other momenta regions which are not satisfied by the model in [11] though 
they argued that this mades only a small numerical difference of the order of 
0.05 X 10~^°. In fact, within the large Nc framework, it has been shown [17] that 
in general for other than two-point functions, to satisfy fully the QCD short- 
distance properties requires the inclusion of an infinite number of narrow states. 

The results in [11] for the Goldstone boson exchanges and for the axial- vector 
exchanges can be found in Table 1 and 3, respectively. 
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Tabic 4: Results quoted in Ref. [11] for the axial- vector exchange depending of 
the /i(1285) and /i(1420) resonances mass mixing. 



Mass Mixing 


10^° X 


No New OPE and Nonet Symmetry 




with M=1.3 GeV 


0.3 


New OPE and Nonet Symmetry 




with M= 1.3 GeV 


0.7 


New OPE and Nonet Symmetry 




with M= Up 


2.8 


New OPE and Ideal Mixing 




with Experimental Masses 


2.2 ± 0.5 



Tabic 5: Full hadronic light-by- light contribution to at 0{Nc). The difference 
between the two results of Refs. [8] and [9] is the contribution of the scalar 
exchange — (0.7± 0.1) x 10~^^. This contribution is not included in Refs. [6,7] 
and [11]. 



Hadronic light-by-light at 0{Nc) 


10^° X 


Nonet Symmetry [6, 7] 
Nonet Symmetry -|- Scalar [8,9] 
Nonet Symmetry [8, 9] 
New OPE and Nonet Symmetry [11] 
New OPE and Ideal Mixing [11] 


9.4 ± 1.6 
10.2 ± 1.9 
10.9 ± 1.9 
12.1 ± 1.0 
13.6 ± 1.5 
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3 Next-to-Leading in the 1/Nc Expansion Con- 
tributions 



At next-to-leading (NLO) in the 1/Nc expansion, the pion loop is the dominant 
one and because the pion mass is not much larger than the muon mass m, one 
expects a contribution of the order of 10~^°. To dress the photon interacting with 
pions, a particular Hidden Gauge Symmetry (HGS) model was used in [6,7] while 
a full VMD was used in [8,9]. The results obtained are -(0.45 ± 0.85) x 10"^° 
in [6] and —(1.9 ±0.5) x 10"^*^ in [8]. Both models satisfy the known constraints 
though start differing at 0{p^) in CHPT. It is also known that the full VMD does 
rather well reproducing higher order terms of CHPT while the special version of 
the HGS used in [6] does not give the correct QCD high energy behavior in some 
two-point functions, in particular it does not fulfill the Weinberg Sum Rules, 
see [8] for more comments. Some studies of the cut-off dependence of the pion 
loop using the full VMD model was done in [8] and showed that their final number 
comes from fairly low energies where the model dependence should be smaller. 
The authors of [11] analyzed the model used in [6, 7] and showed that there 
is a large cancellation between the first three terms of an expansion in powers 
of (m^/Mp)^ and with large higher order corrections when expanded in CHPT 
orders but the same applies to the 7r° exchange as can be seen from Table 6 in 
the first reference in [1] by comparing the WZW column with the others. The 
authors of [11] took (0 ± 1) x 10"-'^'^ as a guess estimate of the total NLO in 
1/Nc contribution. This seems too simply and certainly with underestimated 
uncertainty. 



4 Comparison Between Different Calculations 

The comparison of individual contributions in [6-10, 16] and in [11] has to be 
done with care because they come from different model assumptions to construct 
the full relevant four-point function. In fact, the authors of [16] have shown that 
their constituent quark loop provides the correct asymptotics and in particular 
the new OPE found in [11]. It has more sense to compare results for a|J^' either 
at leading order or at next-to-leading order in the l/N^ expansion. The recent 
results for aj^' at LO in the l/N^. expansion is what is shown in Table 5. The nice 
agreement between them within the quoted uncertainty leads us [1] to take 

as a robust result for the hadronic light-by-light contribution to muon anomaly 
at LO in the 1/Nc expansion. 

The results for the final hadronic light-by-light contribution to quoted 
in [6,7], [8,9] and [11] are in Table 6. The apparent agreement between [6,7] 
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Table 6: Results for the full hadronic light-by-light contribution to a 



Pull Hadronic Light-by-Light 


10^° X 


[6,7] 


8.9± 1.7 


[8,9] 


8.9 ± 3.2 


[11] 


13.6 ± 2.5 



and [8, 9] hides non-negligible differences which numerically almost compensate 
between the quark-loop and charged pion and kaon loops. Notice also that [6, 7] 
didn't include the scalar exchange. Comparing the results of [8,9] and [11], as 
discussed above, we have found several differences of order 1.5 x 10^^*^ which 
are not related to the new short-distance constraint used in [11]. The different 
axial-vector mass mixing accounts for —1.5 x 10~^°, the absence of the scalar 
exchange in [11] accounts for —0.7 x 10~^° and the use of a vanishing NLO in 1/A^c 
contribution in [11] accounts for — 1.9x 10~^°. These model dependent differences 
add up to —4.1 x 10^^° out of the final —5.3 x 10~^° difference between [8,9] 
and [11]. Clearly, the new OPE constraint used in [11] accounts only for a small 
part of the large numerical final difference. 



5 Conclusions 

We observe a nice agreement, see Table 5, between the recent model calculations 
of the hadronic light-by- light contribution to at LO in the 1/Nc expansion, 
hence concluding that 

a;,^''^^ = (11±4) X 10"^° (9) 

is a very solid result. We also understand the origin of the final numerical differ- 
ence between the results quoted in [11] and [8,9]. Its origin is not dominated by 
the new OPE constraint found in [11] and it rather comes from the addition of 
several model dependent differences of order 1.5 x 10~^° as discussed above. 

It is possible and desirable to make a new calculation of aj^' using the tech- 
niques developed in [13,17,18] and the new OPE results [11]. 

The authors of [2] have done a conservative analysis of the present situation 
of the hadronic hght-by-hght contribution to including the NLO in the 1/Nc 
expansion contribution. 

Very valuable information about various pieces of the theoretical models used 
to calculate the hadronic light-by- light contribution to can be obtained by 
measuring the 7r° — >• 77*, 7r° — >• 7*7* and 7r° — > e+e~ decays which constrain the 
off-shell 7r°7*7* and 71^7*7 form factors and the subleading //-dependent terms 
discussed in Section 2 and by measuring the 7*7* —>■ 7r"'"7r~, e'^e~ — > tt'^tt" 
processes which constrain the 7r'*"7r~7*7* vertex which dominates the uncertainty 
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of the pion loop contribution. The 77 programme at the upgraded DA$NE-2 
facihty at Prascati is very well suited for these measurements. 
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